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We study spin dependent transport through a magnetic bilayer graphene nano junction configured 
as two dimensional normal/ferromagnetic/normal structure where the gate- voltage is applied on the 
layers of ferromagnetic graphene. Based on the fourband Hamiltonian, conductance is calculated by 
using Landauer Butikker formula at zero temperature. For parallel configuration of the ferromag- 
netic layers of bilayer graphene, the energy band structure is metallic and spin polarization reaches 
to its maximum value close to the resonant states, while for antiparallel configuration, the nanoj unc- 
tion behaves as a semiconductor and there is no spin filtering. As a result, a huge magnetoresistance 
is achievable by altering the configurations of ferromagnetic graphene especially around the band 
gap. 
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I. INTRODUCTION 

Since spin-orbit coupling in graphene 1 is very weak 2 
and also there is no nuclear spin 3 , spin flip length is so 
long about lfim in dirty samples and room temperature^. 
Clean samples are expected to have longer spin co- 
herency. This is good opportunity for spintronic ap- 
plications based on candidatory of graphene. On the 
other hand, graphene has not intrinsically ferromagnetic 
(FM) properties, however, it is possible to induce ferro- 
magnetism externally by doping and defects^, Coulomb 
interactions^ or by applying an external electric field 
in the transverse direction in nanoribbons 7 . Recently, 
HaugeiP proposed the FM correlations due to strong 
proximity of magnetic states close to graphene. The over- 
lap between the wave functions of the localized magnetic 
states in the magnetic insulator and the itinerant elec- 
trons in graphene induces an exchange field on itiner- 
ant electrons in graphene giving rise spin splitting of the 
transport. The exchange splitting which is induced by 
FM insulator Euo in graphene was estimated to be of 
order of 5meV. This splitting which is effectively simi- 
lar to a Zeeman interaction has so large magnitude that 
can have important effects. Such spin splitting can be 
directly evaluated from the transmission resonances or 
magnetoresistance of FM graphene junction. The ferro- 
magnetism leads to a spin splitting effectively similar to 
a Zeeman interaction but of much larger magnitude. The 
induced exchange field is tunable by an in-plane external 
electric fielcP. The possibility of controlling spin conduc- 
tance in FM monolayer graphene insulator has also been 
studied by YokoyamgM It was found that the spin con- 
ductance has an oscillatory behavior in terms of chemical 
potential and the gate voltage. 

Bilayer graphene on the other hand has shown to have 
interesting properties for application in nanoelectronic 
devices such as transistors based on graphene substrate. 
The new type of integer quantum Hall effects^Mand also 
the electronic band gap controllable by vertically applied 
electric field are of its unusual properties in compared 




FIG. 1: Schematic view of normal/ferromagnetic/normal bi- 
layer graphene. Two gate electrodes can be coated on top of 
the magnetic insulator strips which are located on the upper 
and lower layers. 



to monolayer graphene^^. Moreover, the parabolic 
band structure close to the Dirac points transforms to a 
Mexican-hat like dispersion when an electric field is ap- 
plied on graphene. Optical measurements and theoretical 
predictions propose a 200 meV gap in bilayer graphene. 
This controllable gap makes bilayer graphene as an ap- 
propriate candidate for spintronic devices. An effective 
two-band Hamiltonian can describe the low energy exci- 
tations of a graphene bilayer in the regime of low barrier 
heights 16 . However, the four-band Hamiltonian is known 
to give a better agreement with both ex perime ntal data 
and theoretical tight-binding calculation i 12 \ 17 [ Very re- 
cently, spin splitting of conductance in bilayer graphene 
has been investigated by using an effective two-band 
Hamiltonian emerging from low energy approximation^. 
This approximation is valid when energy of electrons hit- 
ting on the potential barrier is about the barrier height. 
Application of a potential difference between upper and 
lower layers intensifies the failure of this approxima- 
tion. On the other hand, magnetoresistance in bilayer 
graphene has been studied in Ref.^by using 8x8 Hamil- 
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tonian when the induced exchange fields are laid in plane 
of each layer with a rotation in their orientations aginst 
each other. By using Kobo- Greenwood formula, they 
have investigated the dependence of conductivity and 
magnet oresistance on temperature and induced exchange 
field. 

Motivated by these studies, based on the four-band 
Hamiltonian and close to the Dirac points, we study 
spin current through a magnetic barrier creating by use 
of the proximity of a ferromagnetic insulator on bilayer 
graphene. Conductance is calculated by use of Landauer- 
Buttiker formula at zero temperature. The parameters 
of the barrier, energy and angle of incident electrons can 
affect transport through a magnetic barrier classifying in 
the propagating or evanescent modes. The dependnce 
of resonant peaks in transmission on system parameters 
is proposed to follow a resonance condtion. We have 
found in some resonant energies and barrier parameters 
and also around the gapped region that a remarkable spin 
polarization and also magnetoresistance can be achieved. 

Model: We consider a normal/ferromagnetic/normal 
bilayer graphene nano junction. The model we have used 
for a bilayer graphene sandwiched between two ferromag- 
netic insulator is shown schematically in Fig. [I] The 
exchange fields induced by the ferromagnetic insulators 
is supposed to be perpendicular to the graphene plane. 
Therefore, Hamiltonian of the spin up detaches from the 
spin down. Two gate electrodes can be attached to the 
ferromagnetic graphene from the upper and lower layers 
which control the barrier height in each of layers. This 
set-up is different from the systems studied by Refs.LLiM)] 
The exchange field splits this potential depending on the 
spin parallel (+) or antiparallel (— ) to the exchange field. 
So in the ferromagnetic part of bilayer graphene, we have 
y ± = Vo =F A where A and Vo are the exchange field and 
the potential barrier made by the gate voltage, respec- 
tively. So the two spins are scattered from the barriers 
with different heights. It means that energy shift of the 
top of the valance band in the barrier is different for par- 
allel and antiparallel spins to the exchange field. This 
spin splitting causes to shift conductance as a function 
of energy for each spin resulting in magnetoresistance. 
To investigate spin polarization and also magnetoresis- 
tance, we have considered two different configurations so 
that the exchange field inducing by the magnetic insula- 
tors on each layer are parallel or antiparallel with respect 
to each other. The configuration prepared for observa- 
tion of magnetore sista nce differs from the configuration 
considered by Ref! 19 * 20 i The parallel configuration has 
a metallic behavior while the antiparallel configuration 
induces a potential difference between upper and lower 
layers concluding that the system has a semiconductor 
behavior with the band gap of 2 A. 

This paper is organized as follows: we briefly explain 
the formalism which is used for calculating of transmis- 
sion based on the four-band Hamiltonian. Before present- 
ing our results, it is so important to have a short review in 
section III on transport through a barrier deposited on bi- 
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FIG. 2: Schematic view of the potential barrier with height 
Vo and width w and the wave number at normal incidence 
directed in the x-axis for three regions. 



layer graphene and its dependence on the system param- 
eters such as energy of incident quasi-particles and their 
angle hitting into the barrier and also barrier parame- 
ters. The method presented in section II is a detailed 
analysis accompanied with some small corrections on the 
method used by Refill We will present spin polarization 
in the parallel configuration in section IV. Magnetoresis- 
tance and its dependence on energy of incident particles 
and also induced magnetic field will be investigated in 
section V. Finally, the last section concludes our results. 



II. FORMALISM 

In the unit cell of bilayer graphene, we suppose that 
two independent sublattices A and B related to each 
monolayer graphene are connected to each other in the 
Bernal stacking. Close to the Dirac points and in the 
nearest neighbor tight binding approximation, the four 
band Hamiltonian and also its eigenfunction is written 
as the following: 



H 



Trt Vi 

t± V 2 7T+ 
\ 7T V 2 J 



W/ 



(1) 



where 



(Px + ip y ) Vf = -ihvf (d x - k y ) 



and in the above formula k y = k sin where 6 and k are 
the incident angle and wave number of quasi-particles 
hitting on a barrier which is created by applying a bias 
to a metallic strip deposited on bilayer graphene. More- 
over, Vi and V 2 are the gate potentials applied on the 
upper and lower layers of bilayer graphene. Such a gate 
potential can be manipulated by applying a perpendic- 
ular electric field on graphene sheet. Here, the barrier 
is approximated by a square potential of barrier with 
sharp variation. By solving the eigenvalue equation of 



3 



H$> = E$>, the four band spectrum can be concluded as 
the following: 



( £ ') 2 = k 2 + S 2 + ^±t\l4k 2 S 2 /(t') 2 + C-) 2 + k 2 (2) 



where the above parameters are defined as the follow- 
ing: 



(E - V )/hv f = a 



- v , V = 
V 2 )/2hv f , 



t'-- 



- V 2 )/2 

t/hv f . 



(3) 



In the case of S = and k « t' in low energy limit, 
energy spectrum behaves as E — Vo = ±h 2 k 2 /2m where 
m = t/(2v 2 ) is an effective mass. This approximation 
which results in a effective two-band Hamiltonian is valid 
when energy of incident electrons is close to the barrier 
height. On the other word, in the case of zero poten- 
tial difference between two layers, the absolute value of 
(E — Vo) should be much smaller than the interlayer cou- 
pling strength (0.4eV). For S ^ 0, this approximation 
may fail for large potential differences. So one should 
care to choose valid energy and potential ranges. How- 
ever, in spite of RefP^l, in this paper, we use four-band 
Hamiltonian in which the only approximation is the Dirac 
cone. 

If we assume plane wave solution for the Schroedinger 
equation, the wave function in each region with a con- 
stant potential is written as the following matrix prod- 
uct. 



= GM 



(4) 



where matrix elements of the matrices G and M are de- 
fined as the following: 
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FIG. 3: Transmission at normal incidence 6 = as a function 
of energy difference between energy of incident particles and 
the barrier height (E — Vo)/Hvf for potential difference of 
upper layer against lower layer (5 = and lOmeV. 



Here a is the wave number in the x direction and is de- 
fined as: 



a 2 ± = [5 2 + (eT-k 2 y ±vWWTWWW^j} (7) 

In the special case of normal incident angle and zero gate 
potential where k y = (5 = 0, the wave number a+ is real 
in the energy range of < e' < t' and a_ is real if 
the energy of incident particles is in the range of e' < 
0,s' > —t'. In this paper, studied system contains a 
magnetic or electrostatic barrier of potential as shown in 
Figj2j The barrier width is w. The electrostatic potential 
which plays the role of the gate voltage is set to be Vo 
in the barrier part and zero in the first and last regions. 
We suppose that the energy range of incident particles 
is limited to the range of < e' < t' . Consequently 
in the barrier part we have — Vo < e' 2 < t' — Vo. The 
wave numbers behind and in front of the barrier and 
are real while and are imaginary. In the 

(2) (2) 

barrier part, for e' 2 < 0, the wave numbers a\_ J and a_ J 
are imaginary and real, respectively, while for e' 2 > 
they behave vice versa. A schematic view of the barrier 
at normal incidence and wave numbers in each part are 
shown in Fig (2] 

By applying continuity of the wave functions on the 
boundaries of the barrier, one can connect coefficient ma- 
trix of the wave function for the last region As to the 
coefficient matrix for the first region A\. 

Al = NA * (R) 
N = M 1 - 1 (0)G 1 - 1 G 2 M 2 (0)M 2 - 1 (w;)G^ 1 G3M3(w) { ' 

where N is called as transfer matrix. Since and 
are imaginary in the interested energy range, that part of 
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FIG. 4: Contour plot of transmission in plane of the incident 
angle and energy difference of e' 2 — (E—Vo)/hvF accompanied 
with the band structure for potential difference of the upper 
layer against the lower layer to be as a) S = and b) lOmeV. 



wave function which are associated by such wave numbers 
are exponentially a growing or decaying function. So we 

have to set the coefficient of plane wave e ta - )x (c in Eq.J^) 
to be zero for the first region, because this part of wave 
function grows exponentially when x — >• — oo. Therefore, 
coefficients matrix in the first region is supposed to be as 
A\ = [l,r, 0, e p ] T , where the superscript T refers to the 
transpose of a matrix and e g is the coefficient of growing 
evanescent state and r is the coefficient of the reflected 
part of the wave function. In the last region, we have 
to set the coefficient (d in Eq. [4]) of e ~ ia -' >x to be zero 
because this part of the wave function increases exponen- 
tially when x — >> 00. Therefore, the coefficients matrix in 
the last regions is supposed to be as A 3 = [t,0,e<i,0] T 
where t is the coefficient of the transmitted part of wave 
function and is the coefficient of decaying evanescent 
state. In this region, there is no reflected wave. How- 
ever, in equation 8 of Ref. 17 , matrices A\ and A 3 have 
been considered to be completely displaced which leas to 
different results. 

By rearrangement of the transfer matrix elements of 
Eq. |8j the coefficient of transmitted part of wave function 
is derived in terms of transfer matrix elements as the 
following; 



t=[N n -N 13 N 31 /N 33 }- 



(9) 



Since the first and last regions possess similar wave num- 
bers, transmission probability is given as T = \t\ 2 . Be- 
fore presenting our results, in the next section, we will 
shortly review transport properties through a potential 
barrier by using the mentioned formalism. 



III. TRANSMISSION THROUGH A BARRIER 
ON BILAYER GRAPHENE 

The Klein tunneling in monolayer graphene results 
in a complete transmission through a barrier poten- 
tial in normal incident. However, in contrast to mono- 
layer graphene, as a result of chiral symmetry in bilayer 
graphene, transmission is zero for quasiparticles with en- 
ergies lower than the barrier height. In the special case 
of 8 = k y = 0, transmission through a potential barrier 
can be analytically calculated in normal incident and for 
two ranges of energy e' < and e' > 0. 



t 



where 



cos(cr 2 )iu) — iQ sin(cr 2 )iu)] 



(10) 



where in the above formula, parameters are defined as 
s'i = s/hvf,s' 2 = (e — Vo)/hvf. So the real part of the 
wave numbers inside and outside of the barrier part are 
defined as a« = [(s^) 2 + e'J'} 1 ' 2 and = [(s' 2 ) 2 + 
e^t'] 1 / 2 , respectively. The energy of incident particle is 
supposed to be always e[ > 0. So aW is always real. 

For the energy range of E < Vo, the wave number in- 
side the barrier part and consequently Q are imag- 
inary so that = w and Q = iq where k and q are 
real. As a result, transmission tends to zero as a function 
of the system parameters such as w and s' 2 . 



T(0 = 0, 4 < 0) = tt* = [cosh 2 (>w) + q 2 sinh 2 (^)]- 1 

This behavior is the trace of chiral symmetry in bilayer 
graphene. However, if the incident angle is nonzero, 
some resonant peaks appear in the transmission curve 
(see Fig. For the energy range of E > Vo, all param- 
eters such as and Q are real. Thus transmission has 
an oscillatory behavior as a function of e' 2 as the following 
form, 



T(<9 = 0,4 > 0) = [cos 2 (a 



( 2 )„ 



Q 2 sin 2 ^™)]- 1 



(12) 



In the high energies limit E ^> Vq, we have Q — > 1 
and so transmission is complete (T — > 1). By apply- 
ing a vertically electric field in the barrier part, a band 
gap is opened in the band structure of bilayer graphene 
which is proportional to the potential difference between 
potentials of each layers. In this case, chiral symme- 
try is failed and therefore transmission in normal inci- 
dence is nonzero for energies lower than the barrier height 
(E < Vo). Transmission at normal incidence is repre- 
sented in Fig. [3] as a function of e' 2 for S = and lOmeV. 
Application of a vertically electric field causes to emerge 
some resonant tunneling states for energies of E < Vq. In 
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FIG. 5: Contour plot of transmission in plane of the incident 
angle and the barrier width for the barrier height 50meV, the 
energy of incident angle 17meV and for the case of S = 0. 
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the exchange field direction in the parallel configuration, b) 
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this energy range, resonant states originates from inter- 
ference of the incident and scattered waves. For all cases 
such as nonzero incident angles and S ^ 0, the resonant 
peaks are interpreted by the resonance condition relation 
proposed as 



ai/w = nn 



(13) 



where a& is the x-component of the wave number inside 
the barrier which can be calculated by Eq. [7J 

To have more complete view, we prepare a contour 
plot of transmission in plane of the incident angle and e' 2 
which is shown in Fig. [4] for a fixed width of the barrier. 
For the normal incidence {6 = 0), transmission behavior 
is compatible with the results shown in Fig. [3j 

For energies higher than the barrier height e' 2 > 0, 
transmitting channels are opened over all ranges of en- 
ergies. However, transmitting window for the incident 



FIG. 7: a) Conductance and b) spin polarization as a func- 
tion of e' 2 for different induced exchange fields A in the par- 
allel configuration. Here, the barrier height and width are 
considered to be as 50 meV and 40 nm, respectively. 



angles is limited with the condition that (in Eq. ^ 
is real. In the case of 5 = 0, the range of incident an- 
gle in which transmission is high can be extracted as 
-1 y/(4) 2 +4^ 



sin 



< < sin 



k 



Therefore, 



by increasing s 2 , the range of angles with high trans- 
mission becomes more extended. In the energy range of 
e 2 < 0, independent of the value of 5, resonant peaks 
emerge for nonzero incident angles (6 ^ 0) which obey 
the resonance condition a^w = utt. So additional to some 
resonant energy states, we have some resonant widths in 
which transmission is high. Fig. [5] shows transmission 
in plane of the incident angle and the barrier width for 
e' < and 5 = 0. It is shown that based on the reso- 



nance condition (Eqs. 13|7 ), in large incident angles, a& 
reduces and so in a fixed resonance order (n), the reso- 
nance condition is satisfied for wide barriers. Therefore, 
the resonance strips with complete transmission shown 
in Fig. [5j depend strongly on the incident angle in the 
range of wide barriers. 

By applying a vertically electric field in the barrier 
part, a band gap is opened around e 2 = 0. This band 
gap also has a trace in transmission as a transport gap 
shown in Fig. Hb. 



IV. RESULTS 

By application of an averaged gate voltage Vb, band 
structure in the barrier part is shifted by Vq value. Fig. [6] 
shows band structure of parallel and antiparallel configu- 
ration magnetic insulators when a gate voltage is applied 
on the barrier part. In case the exchange fields inducing 
in each layers of bilayer graphene are parallel, particles 
with spin parallel (spin up) and antiparallel (spin down) 
to the exchange fields are scattered from barriers with dif- 
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FIG. 8: a) Conductance and b) spin polarization as a func- 
tion of barrier width for different induced exchange fields A 
in the parallel configuration. Here, the barrier height and in- 
cident energy are considered to be as 50 meV and 17 meV, 
respectively. 



ferent heights. In the parallel configuration, spin splitting 
of the barrier potential in the ferromagnetic graphene is 
written as V~ — V + = 2 A. Such spin splitting is also seen 
in the band structure that is shown in Fig. |6^l. It is seen 
that the top of valance band are shifted to lower/higher 
energies for spins up/down. However, in the antiparallel 
configuration, the band structure shown in Fig. ^jp is the 
same for both up and down spins. A band gap which is 
proportional to 2A appears in the band structure of the 
antiparallel configuration. 



A. Spin Polarization 

Here, there is a correspondence between the band 
structure and transmission. According to the band struc- 
ture, we expect to emerge spin polarization just for par- 
allel configuration because energy bands for up and down 
spins are shifted by 2 A with respect to each other. How- 
ever, since the band structure for antiparallel configu- 
ration is the same for both spins, it is not expected to 
have spin polarization for this configuration. The spin 
polarization is defined as: 



Gdo 



Gup H~ Gdown 



(14) 



where G up and Gdown are conductance for up and down 
spins. The conductance is calculated by using Landauer 
formalism in the linear regime. Therefore, conductance 
is proportional to angularly averaged transmission pro- 
jected along the current direction. 



pir/2 

G= T(E,cos(0))cos0d0 

J—k/2 



It is clear that additional to the transmission curves 
(FigQ, resonance peaks also appears in conductance. 
Since up and down spins in the parallel configuration see 
barriers with different heights, resonance peaks in con- 
ductance as a function of Fermi energy E are shifted to 
higher energies as A for spin down and to lower energies 
as —A for spin up. This mismatching of conductance 
peaks for two spins causes to a large spin polarization at 
resonance states. Fig. [7] displays conductance and spin 
polarization as a function of e' 2 for the parallel configu- 
ration. It is shown that conductance peaks and conse- 
quently spin polarization appears in the energy range of 
e 2 < 0. It is seen that by inducing an exchange field, con- 
ductance peaks in Fig. [7^ split into two peaks which are 
related to each spin. This spin splitting is about 2 A. Spin 
polarization shown in Fig. [7}} has an oscillatory behavior 
with energy of incident particles for energies lower than 
the barrier height e 2 < 0. The amplitude of spin polar- 
ization increases with the induced exchange field A and 
reaches to its maximum value. However, spin polariza- 
tion tends to zero for energies greater than the potential 
height e 2 > except at E ~ Vq. 

In the parallel configuration and for e' 2 < 0, Fig. [8^l 
shows that conductance in the resonance widths has a 
peak. These peaks which are also seen in the transmission 
curves of Fig. [5] are explained by the resonance condition 
of Eq. 13 It is shown that spin splitting of conductance 



peaks also appears in the resonance widths which is orig- 
inated from different barrier heights for two spins up and 
down. It should be noted that the conductance at res- 
onance widths decreases for wide barriers. In the wide 
range of widths, the angularly window for transmitting 
channels shown in Fig. [5] decreases with the widths. 

Fig. [8}} shows spin polarization as a function of the 
barrier width. Again, spin polarization has an oscilla- 
tory behavior with the barrier width. The amplitude of 
spin polarization strongly increases by an increase of the 
induced exchange field. Therefore, to manifest this spin 
polarization, we should manufacture the ferromagnetic 
graphene part with the special widths in which spin po- 
larization reaches to the value of unity. 



B. Magnetoresistance 

In this section, we will show that by switching between 
parallel and antiparallel configurations, one can obtain 
large magnetoresistance. Magnetoresistance is defined as 
the following: 



MR 



qv _ c a P 

Qp + Qap 



(15) 



where G p = G% + G p down and G ap = G% + Gf own are 
conductance for parallel and antiparallel configurations. 

Fig. [9] displays conductance in the parallel and antipar- 
allel configurations and also magnetoresistance as a func- 
tion of e' 2 and the barrier width for a fixed exchange field 
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FIG. 9: Conductance in the parallel and antiparallel config- 
urations as a function of a)^ = (E — Vo)/Hvf for a barrier 
with the width of 40 nm and, c) the barrier width for a barrier 
with the height of 50 meV. Magnetoresistance as a function 
of b) e' 2 for a barrier with the width of 40 nm, d) the barrier 
width for a barrier with the height of 50 meV. The induced 
exchange field is considered to be as A = 5meV. 



A = 5meV. As we before expressed, spin splitting at 
the resonance states (for e' 2 < 0) emerges in conductance 
peaks in the case of the parallel configuration. This be- 
havior is clear in Fig. [9^i and[9]3. However, this splitting 
will not occur for the case of antiparallel configuration. 
Therefore, large magnetoresistance appears around the 
conductance resonance peaks. In the parallel configura- 
tion, a band gap appears around the barrier edge in the 
interval Vq — A < E < Vo + A. This band gap has a trace 
in transmission and consequently conductance. Zero con- 
ductance region around the barrier edge e 2 ~ which is 
seen in Figj9^i, is a result of the band gap. Since there is 
no such a band gap in the parallel configuration, magne- 
toresistance as shown in Figj9}} reaches to its maximum 
value in the energy band gap. In the energy range greater 
than the barrier height e' 2 > 0, there is no spin splitting 
and therefore, magnetoresistance tends to zero. 

As we showed before, conductance has peak at reso- 
nant widths. Similar to the previous case, spin splitting 
occurs just for the parallel configuration. So magnetore- 
sistance increases around the resonance widths. The os- 



FIG. 10: a) Conductance in the parallel and antiparallel con- 
figuration and b) magnetoresistance as a function of the in- 
duced exchange field A for a barrier with the height of 50 
meV and energy of incident particles as 40 meV. Here the 
barrier width is 20 nm. 



cillatory behavior of magnetoresistance as a function of 
the barrier width is represented in Figj9]i. 

As we showed, there is a large magnetoresistance 
around the barrier edge E w Vq. In this range of ener- 
gies, we investigate the dependence of magnetoresistance 
to the induced exchange field. This exchange field of 
graphene can be controlled by an in-plane external elec- 
tric fielcP. Fig. 10 d represents that magnetoresistance 
increases monotonically by increasing the exchange field 
A. It is interesting by increasing the exchange field up 
to lOmeF, magnetoresistance reaches to its maximum 
value. 



To explain this behavior, we investigate the depen- 
dence of conductance on the exchange field in the par- 
allel and antiparallel configurations. In the antiparal- 
lel configuration, the band gap which is limited in the 
interval of Vb — A < E < Vq + A, enhances by in- 
creasing the exchange field. Therefore, conductance in 
the antiparallel configuration goes to zero when the ex- 
change field is increased. Suppression of the conductance 
with the exchange field in the antiparallel configuration 
is shown in Fig. 10 1. However, in the parallel config- 
uration, conductance increases by enhancement of the 
exchange field. The reason of this enhancement comes 
back to have larger angularly transmitting windows for 
larger e 2 (see Fig.[4|. In fact, effective potential for spins 
up V + = Vq — A is decreased by an increase in A. So 
e' 2 = (E — Vo)/hvp for a fixed energy is increased and 
consequently G up and so G p is increased by A. As a 
conclusion, for the exchange fields up to 10 meV, sup- 
pression of G ap and an increase of G p results in a large 
magnetoresistance which is so useful for designing spin 
memory devices. 



8 



V. CONCLUSION 

We have studied spin polarization and magnetoresis- 
tance of a normal/ferromagnetic/normal junction of bi- 
layer graphene by using transfer matrix method and 
based on the four-band Hamiltonian. Transport prop- 
erties simultaneously is controlled by two gate electrodes 
(Vo), which are applied on the ferromagnetic graphene. 
Two configurations of the exchange field is considered 
perpendicular to the graphene sheet. This exchange field 
is induced by the proximity of a localized magnetic or- 
bital in a magnetic insulator coating on top of each layers 
of bilayer graphene in the barrier part. In the parallel 
configuration which graphene has a metallic behavior, a 
spin splitting 2A occurs for the conductance at the reso- 



nant states just for energies lower than the barrier height 
E < Vo- However, there is no spin splitting in the an- 
tiparallel configuration. A band gap of 2A is opened 
in the antiparallel configuration which makes it a semi- 
conductor. As a result of spin splitting in the parallel 
configuration, an oscillating spin polarization emerges for 
energies lower than the barrier height. Furthermore, an 
oscillatory of magnetoresistance with large amplitude is 
achievable for E < Vq when we are able to switch between 
two configurations. There is also a large magnetoresis- 
tance in the energy range around the barrier edge origi- 
nating from the band gap which is openned by a vertically 
electric field. In this range of energy, magnetoresistance 
reaches to its maximum value when the exchange field is 
increased by an in-plane external electric field. 
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